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Abstract 

Krawtchouk's polynomials occur classically as orthogonal polyno- 
mials with respect to the binomial distribution. They may be also 
expressed in the form of matrices, that emerge as arrays of the values 
that the polynomials take. The algebraic properties of these matri- 
ces provide a very interesting and accessible example in the approach 
to probability theory known as quantum probability. First it is noted 
how the Krawtchouk matrices are connected to the classical symmet- 
ric Bernoulli random walk. And we show how to derive Krawtchouk 
matrices in the quantum probability context via tensor powers of the 
elementary Hadamard matrix. Then connections with the classical 
situation are shown by calculating expectation values in the quantum 
case. 



1 Introduction 



Some very basic algebraic rules can be expressed using matrices. Take, for 
example, 



(a + 6)2 

(a + b){a - h) 

(a - hf 



a? + 2ah + 
- 2ab + 



1 1 

2 
1 -1 



(the expansion coefficients make up the columns of the matrix). In general, 
we make the definition: 

Definition 1.1 The A^*'^-order Krawtchouk matrix K^^'> is an (A^ + 1) x 
(A^ + 1) matrix, the entries of which are determined by the expansion: 



N 



(N) 



(1) 



1=0 



The left-hand-side G{v) = (l + v)^^^ (1 —vY is thus the generating function 
for the row entries of the j^^ column of K^^\ Expanding gives an explicit 
expression: 



K^P - E(-i) 



Ms 



k 



Here are the Krawtchouk matrices of order zero and one: 



i^(o) 



1 1 
1 -1 



(2) 



More examples can be found in Table 1 of Appendix 1. In the remaining of 
the text, matrix indices run from to A'^. 



One may view the columns of Krawtchouk matrices as generalized bino- 
mial coefficients. The rows define Krawtchouk polynomials: for a fixed order 
A?", the i-th Krawtchouk polynomial is the function 

K,{j,N)^k\P 

that takes its corresponding values from the i-th row. One can easily show 
that Ki{j, N) is indeed a polynomial of degree i in the variable j. 
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Historically, Krawtchouk's polynomials were introduced and studied by 
Mikhail Krawtchouk in the late 20's [T^llIS]- Since then, they have appeared 
in many areas of mathematics and applications. As orthogonal polynomials, 
they occur in the classic work by Szego ^2^. They have been studied from 
the point of view of harmonic analysis and special functions, e.g., in work 
of Dunkl [HI El- In statistical considerations, they arose in work of Eagleson 
[8j and later Vere- Jones [21]. They play various roles in coding theory and 
combinatorics, for example, in Mac Williams' theorem on weight enumerators 
[T71 [14j . and in association schemes [31 HI [5]. 

A classical probabilistic interpretation has been given in [TU] . In the con- 
text of the classical symmetric random walk, it is recognized that Krawtchouk's 
polynomials are elementary symmetric functions in variables taking values 
±1. Specifically, if C,i are independent Bernoulli random variables taking 
values ±1 with probability |, then if j of the are equal to —1, the i^^ 

elementary symmetric function in the is equal to k[^\ It turns out that 
the generating function ([1]) is a martingale in the parameter N . Details are 
in Section [3] below. 

As matrices, they appeared in the 1985 work of N. Bose [IJ on digital 
filtering, in the context of the Cayley transform on the complex plane. The 
symmetric version of the Krawtchouk matrices has been considered in [9] . 

Despite this wide research, the full potential, meaning and significance 
of Krawtchouk polynomials is far from being complete. In this paper we 
look at Krawtchouk matrices as operators and propose two new ways in 
which Krawtchouk matrices arise: via classical and quantum random walks. 
Especially the latter is of current interest. The starting idea is to represent 
the second Krawtchouk matrix (coinciding with the basic Hadamard matrix) 
as a sum of two operators 



" 1 1 ' 




"or 




" 1 


" 


1 -1 




1 


+ 





-1 



Via the technique of tensor products of underlying spaces we obtain a rela- 
tionship between Krawtchouk matrices and Sylvester-Hadamard matrices. 

The reader should consult Parthasarathy's [18] for material on quantum 
probability. It contains the operator theory needed for the subject as well as 
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showing the connections with classical probability theory. 

For information on Hadamard matrices as they appear here, we recom- 
mend Yarlagadda and Hershey's work [22] which provides an overview of 
the subject of Sylvester-Hadamard matrices, indicating many interesting ap- 
plications. For statisticians, they point out that in Yates' factorial analysis, 
the Hadamard transform provides a useful nonparametric test for association. 

Yet another area of significance of this research lies in the quantum com- 
puting program [151 [H] • Details on this connection will appear in an inde- 
pendent work. 

This paper is organized as follows. In the next section, we review basic 
properties of Krawtchouk matrices. The identities presented, although basic, 
seem to be new and do not appear in the references cited. Section [3] presents 
the classical probability interpretation. It may be viewed as a warm-up lead- 
ing to the quantum random walk introduced and studied in Section HI and 
to the relationship between Krawtchouk matrices and Sylvester-Hadamard 
matrices. The generating function techniques used there are original with 
the present authors. In the last subsection, calculating expectation values 
in the quantum case shows how the quantum model is related to the classi- 
cal random walk. Appendix 1 has examples of Krawtchouk and symmetric 
Krawtchouk matrices so that the reader may see concretely the subject (s) 
of our discussion. Appendices 2 (tensor products) and 3 (symmetric tensor 
spaces) are included to aid the reader as well as to clarify the notation. 



2 Basic properties of Krawtchouk matrices 

(1) The square of a Krawtchouk matrix is proportional to the identity matrix. 

(ir(^))2 = 2^/. 

This remarkable property allows one to define a Fourier-like Krawtchouk 
transform on integer vectors. 

(2) The top row is all I's. The bottom row has il's with alternating 
signs, starting with +1. The leftmost entries are just binomial coefficients. 
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K^^"^ — . The rightmost entries are binomial coefficients with alternat- 



jAN) 

mg signs, Kl^' 



(3) There is a four-fold symmetry: \K^ 



1 3 I 



\^N-ij\ 



\k: 



(N) 



i N-j\ 



N-j\ 



Krawtchouk matrices generalize Pascal's triangle in the following sense: 
Visualize a stack of Krawtchouk matrices, the order N increasing downwards. 
Pascal's triangle is formed by the leftmost columns. It turns out that Pascal's 
identity holds for the other columns as well. Less obvious is another identity 
— call it dual Pascal. 



Proposition 2.1 Set a = i^S^-, b = Kfy , A = Kf]^", B = Kf];^ . 

1. (Cross identities) The following mutually inverse relations (Pascal 
and dual Pascal) hold: 



(N) 



a + b^ A 
b-a^B 



and 



A + B^2b 
A-B^2a. 



2. (Squcire identity) In a square of any four adjacent entries in a 
Krawtchouk matrix, the entry in the left-bottom corner is the sum of the 
other three, i.e.. 



for K — 



a c 
b d 



one has b — a + c-\- d. 



Proof. For a + 6, consider 

(1 + vf+'-^{l - vy = (1 + (1 + v)'^~^{l - vy . 
For b — a, consider 

(1 + - vy+^ = (1 - (1 + - vy . 

The inverse relations are immediate. The square identity follows from the 
observation {a + b) + {c + d) — A + B — 2b, hence a + c + d — b. u 
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The square identity is useful in producing the entries of a Krawtchouk 
matrix: fiU the top row with I's, the right-most column with sign- alternating 
binomial coefficients. Then, apply the square identity to reproduce the ma- 
trix. 

In summary, the identities considered above can be written as follows: 



Cross identities: 

(0 -H-lj 



(JV) 



.(iV-l) 



(N) 



Squaire identity: 



(TV) _ r^(N) r^(N) 



j+1 ^ 



If each column of the matrix is multiplied by the corresponding binomial 
coefficient, the matrix becomes symmetric. Let B^^^ denote the (A'" -|- 1) x 
[N + 1) diagonal matrix with binomial coefficients 



. i J 



(3) 



as its non-zero entries. Then, for each N > 0, one defines the symmetric 
Krawtchouk matrix as 



Example: For A?" = 3, we have 



1111 
3 1-1-3 
3-1-1 3 
1-1 1-1 



10 

3 

3 

1 



1 3 
3 3 
3 -3 
1 -3 



3 1 

-3 -3 

-3 3 

3 -1 



Some symmetric Krawtchouk matrices are displayed in Table 2 of Appendix 
1. 
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3 Krawtchouk matrices and classical random 
walk 

In this section we will give a probabilistic meaning to the Krawtchouk ma- 
trices and some of their properties. 

Let C,i be independent symmetric Bernoulli random variables taking values 
±1. Let X^v = ^1 + ■ ■ ■ + ^Af be the associated random walk starting from 
0. Now observe that the generating function of the elementary symmetric 
functions in the C,i is a martingale, in fact a discrete exponential martingale: 

N 

= n(i + = E ^'"4^1' • • • , Civ) , 

i=l k 

where ak denotes the k^^ elementary symmetric function. The martingale 
property is immediate since each has mean 0. Suppose that at time N, 
the number of the that are equal to —1 is Jn, with the rest equal to +1. 
Then = {N — Xn)/2 and can be expressed solely in terms of N and 
Xjv, or, equivalent ly, of N and Jn 

Mn = {1 + Vf-^^{1 - VY^ = (1 + ^)(^+^iv)/2(l _ ^)(A'-^iv)/2 _ 

From the generating function for the Krawtchouk matrices, ([1]), follows 



'■Jn ' 



SO that as functions on the Bernoulli space, each sequence of random vari- 
ables K^j^ is a martingale. 

Now we can interpret two basic recurrences of Proposition l2.1[ For a fixed 
column of K^^\ the corresponding column in satisfies the Pascal 

triangle recurrence: 

To see this in the probabilistic setting, write Mat+i = {1 + v^n)Mn . Observe 
that for to remain constant, must take the value and expanding 
(1 + v)M]^ yields the Pascal recurrence as in the proof of Proposition 12.11 It 
is interesting how the martingale property comes into play. We have 



(N+l) 
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since half the time ^at+i is —1, increasing by 1, and half the time is 
unchanged. Thus, writing j for jjv, 

Af' = I ««" + «>) . 

Many further properties of Krawtchouk polynomials may be derived from 
their interpretation as elementary symmetric functions on the Bernoulli space 
with scope for probabilistic methods as well. 



4 Krawtchouk matrices and quantum random 
walk 

In quantum probability, random variables are modeled by self-adjoint opera- 
tors on Hilbert spaces and independence by tensor products. We can model 
a symmetric Bernoulli random walk as follows. Consider a 2-dimensional 
Hilbert space ^ = and two special 2x2 operators. 



1 

1 



and G 



1 
-1 



satisfying F'^ — G'^ — I (the 2x2 identity). The fundamental Hadamard 
matrix H coincides with the second Krawtchouk matrix. Now we shall view 
it as a sum of the above operators 



H^F + G 



One can readily check that 

FH = F{F + G) = {F + G)G = HG 



(4) 



(use F^ = G^ = I). This, of course, can be viewed as the spectral decompo- 
sition of F and we can interpret the Hadamard matrix as a matrix reducing 
F to diagonal form. 

Remcirk 4.1 Note that the exponentiated operator 

( p^ _ cosh^ sinh^ 
^"^^ sinhz coshz 
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has the expectation value in the state Cq equal to 

(eo, exp(2F)eo) = coshz, (5) 

where Cq denotes the transpose of [1, 0]. This coincides with the moment gen- 
erating function for the symmetric Bernoulli random variable taking values 
±1, showing that indeed we are dealing with the (quantum) generalization 
of the classical model. 



The Hilbert space of states is represented by the A^-th tensor product of 
the original space V, that is, by the 2^-dimensional Hilbert space V^®^ (see 
Appendix 2 for notation). Define the following linear operators, in all, in 

/2 = J®F®/®---®/ 



each fi describing a flip at the z-th position. These are the quantum equiv- 
alents of the random walk variables from Section [3l We shall consider the 
superposition of these independent actions, setting 



= /i + ■ ■ ■ + /, 



N 



Notation: For notational clarity, since is flxed throughout the discussion, 
we drop the index N from the X's. 

Analogously, we define: 

g2 = I (g) G ® I ■ ■ ■ ® I 

qn = /®/®---®G, 

with Xg = gi + ■ ■ ■ + Qn- Finally, let us extend H to the A^-fold tensor 
product, setting iJ^r = H®'^ . These are the well-known Sylvester-Hadamard 
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matrices with the first few given here: 



• o«o«o«o 

• •oo««oo 

• oo«*oo« 

• •••oooo ' 

• o«oo«o« 

• •oooo«« 

• oo«o««o 

etc., where, for typographical reasons, we use • for 1 and o for —1. 

It turns out that our X-operators intertwine the Sylvester-Hadamard 
matrices. For illustration, consider a calculation for N = 3: 

where the relation FH — HG is used. This clearly generahzes to /fei^Ar = 
Hisigk and, by summing over k, yields an important relation: 

XfHn = H^Xq . 

Now, we shall consider the symmetrized versions of the operators (the reader 
is referred to Appendix 3 for the theory and methods used here). Since 
products are preserved in the process of reduction to the symmetric tensor 
space, we get 

XfHn — HnXq , 

the bars indicating the corresponding induced maps. We know how to cal- 
culate H N from the action of H on polynomials in degree N . For symmetric 
tensors the components in degree arc 

where < ^ < A^. 

Proposition 4.2 For each N > 0, symmetric reduction of the N^^ Hadamard 
matrix results in the transposed N^^ Krawtchouk matrix: 

{HN)ij = K^jP . 



Hi 



Ho 



• 


• 


• 


• 


• 


o 


• 


o 


• 


• 


o 


o 


• 


o 


o 


• 



H. 
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Proof. Writing {x,y) for (xo,Xi), we have in degree N for the k^^ compo- 
nent: 



{x + y) 



N-k, 



X 



-yf = Y.Hkix''-W 



Scahng out x^ and replacing v = y/x yields the generating function for the 
Krawtchouk matrices with the coefficient of equal to Kil^\ Thus the re- 
sult. ■ 



Now consider the generating function for the elementary symmetric func- 
tions in the quantum variables fj. This is the A'"-fold tensor power 



T^{t) = {I + tFf^ = I'^^ + tXF + 

noting that the coefficient of t is Xp. Similarly, define 

gj^(t) = (/ + tG)®^ = 7®^ + iXG + --- . 

From (/ + tF)H = H{I + tG) we have 

J'nHn — HnQn and TnHn — H^Qn . 

The difficulty is to calculate the action on the symmetric tensors for opera- 
tors, such as Xp, that are not pure tensor powers. However, from J-N{t) and 
^jv(i) we can recover Xp and Xq via 



X. 



dt 



t=o 



(/ + tF) 



d 
di 



t=o 



(/ + toy 



with corresponding relations for the barred operators. Calculating on poly- 
nomials yields the desired results as follows. 



I + tF 



1 t 
t 1 



I + tG^ 



1 + t 
1-t 



In degree A'^, using x and y as variables, we get the k*^ component for Xp 
and Xg via 



d 
di 



t=o 



(x + tyr-'{tx + yY 



and since I -\-tG is diagonal, 

d 

dt t=o 



{N - 2k) x^-^y^ 
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For example, calculations for = 4 result in 






4 








" 




■ 1 


4 


6 


4 


1 


1 





3 










1 


2 





-2 


-1 





2 





2 





, Hi = 


1 





-2 





1 








3 





1 




1 


-2 





2 


-1 











4 







1 


-4 


6 


-4 


1 



4 

















2 



































-2 

















-4 



Since Xq is the result of diagonalizing Xp, we observe that 

Corollary 4.3 The spectrum of Xp is N, N — 2, ... ,2 — N, —N, coinciding 
with the support of the classical random walk. 



4.1 Expectation values 

To find the probability distributions associated to our Xp operators, we must 
calculate expectation values, cf. Remark [4. II In the present context, expec- 
tation values in two particular states are especially interesting. Namely, in 
the state Cq and in the normalized trace, which is the uniform distribution on 
the spectrum. In the A^-fold tensor product, we want to consider expectation 
values in the ground state | 000 . . . ) and normalized traces. Then we can 
go to the symmetric tensors. 

The scalar product on the tensor product space factors, correspond- 
ing to independence in classical probability. Thus, from ([5]) one obtains 
the expectation value of exp{zXF) in the ground state |000...0) to be 
(coshz)^. Similarly, the trace of the tensor product of operators is the 
product of their traces. So, for the trace, tr exp(zF) = 2 cosh z implies 
tr exp{zXp) = 2^(coshz)^ and, after normalizing, this yields (coshz)^. 
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For the barred operators, we consider the symmetric trace. Here we use 
the symmetric trace theorem, detailed in Appendix 3. It tells us that 
the generating function for the symmetric traces of any operator A in the 
various degrees is 

det(/ — tA)"^. Taking A = exp(zF), we have 

det(/-te^^)-^ = [(l-te^)(l-te-^)]-^ 
= {1 -2t cosh z + t'^y^ . 

The latter is the generating function for Chebyshev polynomials of the second 
kind, Un, so that the normalized symmetric trace is 

{N + iyhrg^exp{zF) = t/^(cosh^)/(iV + 1) , 

which equals as well 

e^(A^+i) - e-^(^+i) _ sinh(iV + l)z 
{e^ - e-^){N + 1) ~ (N+l) sinhz ' 

This corresponds to a uniform distribution on the support of the random 
walk at time N, namely, -TV, 2 - N, . . . , N - 2, N. 

Acknowledgment. We would hke to thank Marios Viana for inviting 
us to participate in the special session and we extend our appreciation for all 
the hard work involved in organizing the session as well as related activities. 
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Appendix 1: Krawtchouk matrices 



i^(3) 



1 1 
1 -1 



1 

1 


1 
1 


1 
1 










9 


n 


_9 










1 


-1 


1 












1 

1 


1 
1 


1 " 
1 








Q 
O 


1 

± 


_1 
— ± 


_Q 








Q 
O 


1 

— ± 


1 

— ± 


Q 
O 








1 
± 


_1 

— ± 


1 
± 


_1 








1 


1 


1 


1 


1 " 






4 


2 





-2 


-4 






6 





-2 





6 






4 


-2 





2 


-4 






1 


-1 


1 


-1 


1 






1 


1 


1 


1 


1 


1 ■ 




5 


3 


1 


-1 


-3 


-5 




10 


2 


-2 


-2 


2 


10 




10 


-2 


-2 


2 


2 


-10 




5 


-3 


1 


1 


-3 


5 




1 


-1 


1 


-1 


1 


-1 




1 


1 


1 


1 


1 


1 


1 


6 


4 


2 





-2 


-4 


-6 


15 


5 


-1 


-3 


-1 


5 


15 


20 





-4 





4 


- 


-20 


15 


-5 


-1 


3 


-1 


-5 


15 


6 


-4 


2 





-2 


4 


-6 


1 


-1 


1 


-1 


1 


-1 


1 



Table 1: Krawtchouk matrices 
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^(4) 



^(6) 



1 2 1 

2 0-2 
1-2 1 



1 


3 


3 1 " 










3 


3 - 


-3 -3 










3 


-3 - 


-3 3 










1 


_3 


3 -1 










1 


4 


6 4 


1 ■ 








4 


8 


-8 


-4 








6 


- 


-12 


6 








4 


-8 


8 


-4 








1 


-4 


6 -4 


1 








1 


5 


10 


10 


5 


1 " 




5 


15 


10 - 


-10 - 


-15 


-5 




10 


10 


-20 - 


-20 


10 


10 




10 


-10 


-20 


20 


10 


-10 




5 


-15 


10 


10 - 


-15 


5 




1 


-5 


10 - 


-10 


5 


-1 




1 


6 


15 


20 


15 


6 


1 


6 


24 


30 


- 


-30 


-24 


-6 


15 


30 


-15 - 


-60 - 


-15 


30 


15 


20 





-60 





60 





-20 


15 


-30 


-15 


60 - 


-15 


-30 


15 


6 


-24 


30 


- 


-30 


24 


-6 


1 


-6 


15 - 


-20 


15 


-6 


1 



Table 2: Symmetric Krawtchouk matrices 
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Appendix 2: Tensor products 



Fulton and Harris ^Jj is a useful reference for this section and the next. Also 
Parthasarathy [18], Chapter II, is an excellent reference. 

Let y be a dimensional vector space over R. We fix an orthonormal 
basis {eo, . . . , es} with d = 1 + 6. Denote tensor powers of V by V"*^^, so that 
Y^2 — Y g|;(3_ j\ basis for V^'^ is given by all A^-fold tensor products of 
the basis vectors Cj , 

I nin2 . . . riAT ) = (g) e„2 (g) ■ ■ ■ (g) e„^ . 

Note that we can label these d'^ basis elements by all numbers to c?^ — 1 and 
recover the tensor products by expressing these numbers in base d, putting 
leading zeros so that all extended labels are of length N. 

Now let {Ai : 1 < i < A^} be a set of linear operators on V. On 
V'^^ , the linear operator A = Ai ® A2 <^ ■ ■ ■ ® A^^ acts on a basis vector 
I nina . . . ) by 

A\ nins . . . ) = AiCn^ (g ■ • ■ (g A^Cnj^ . 
This needs to be expanded and terms regrouped using bilinearity. 

If A and B are two d x d matrices, the matrix corresponding to the 
operator A^ B is the Kronecker product, a d'^ x d"^ matrix having the block 
form: 

aoo-B • • • cbosB 
Oio-B . . . aisB 

asoB . . . assB 

This, iteratively, is valid for higher-order tensor products (associating from 
the left by convention). The rows and columns of the matrix of a linear 
operator acting on V®'^ are conveniently labeled by associating to each basic 

N 

tensor | nin2 ■ ■ ■ ) the corresponding integer label n^d^'^, which thus 

k=l 

provides a canonical ordering. 



16 



Appendix 3: Symmetric tensor spaces 



Here we review symmetric tensor spaces as spaces of polynomials in com- 
muting variables. This material is presented with a view to the infinite- 
dimensional case in p3], pp. 105ff., however we focus on the finite-dimensional 
context and include as well an important observation contained in Theorem 

1131 

The space V"*^^ can be mapped onto the space of symmetric tensors, V"®^^ 
by identifying basis vectors (in \/®^) that are equivalent under all permu- 
tations. Alternatively, one can identify the basic tensor | nin2 . . . ) with 
the monomial x„^x„2 ■ ■ -Xn^ in the commuting variables Xq, . . . ,xs. Hence 
we have a linear map from tensor space into the space of polynomials, itself 
isomorphic to the space of symmetric tensors: 

N>0 N>0 

In the symmetric tensor space, tensor labels need to count only occupancy, 
that is, the number of times a basis vector of V occurs in a given basic tensor 
of V"®^. We indicate occupancy by a multi-index which is the exponent of 
the corresponding monomial. The dimension of V^^'^ is thus 

that is, the number of monomials homogeneous of degree N. 

Given an operator A on V, let = A®'^ . Then An induces an oper- 
ator An on V^^^ from the action of A on polynomials, which we call the 
symmetric representation of A in degree N. For convenience we work 
dually with the tensor components rather with the action on the basis vec- 
tors. Denote the matrix elements of the action of An by Amn- If A has 
matrix entries Aij, let 

Vi = ^ ^ijXj ■ 
3 

Then the matrix elements of the symmetric representation are defined by the 
relation (expansion): 

„»"0 . . . —\^~A ™no . . . rr'^S 
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with multi-indices m and n. 



Composition of Ai with A2 shows that mapping to the symmetric repre- 
sentation is an algebra homomorphism, i.e., 

A1A2 = A1A2 . 

Explicitly, in basis notation 

r 

Define the symmetric trace in degree iV of A as the trace of the matrix 
elements of A^, i.e., the sum of the diagonal matrix elements: 

\m\=N 

with \m\ denoting, as usual, the sum of the components of m. Observe 
that if A is upper-triangular, with eigenvalues Ai, . . . , A^, then the trace of 
this action on the space of polynomials homogeneous of degree N is exactly 
hfyf{Xi, . . . , Arf), the A^*^ homogeneous symmetric function in the A's. 



We recall a useful theorem on calculating the symmetric trace. Since the 
mapping from A to A]\f is a homomorphism, a similarity transformation on 
A extends to one on Ajsf thus preserving traces. Now, any matrix is similar 
to an upper-triangular one with the same eigenvalues, thus follows [19]: 

Theorem 4.4 Symmetric trace theorem Denoting by ti^^^ the trace of 
the symmetric representation on polynomials homogeneous of degree N, 

-< 00 

E ^ Sym^ • 



det{I-tA) 

Proof. With {Aj} denoting the eigenvalues of A, 



1 _ 1 



det{I -tA) 



N=0 



as stated above. 



Sym^ 
N=0 
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Remark 4.5 Note that this result is equivalent to MacMahon's Master The- 
orem in combinatorics [TB] . 

Remark 4.6 From another point of view, Chen and Louck [2j, considering 
powers of the bilinear form J^ij ^i^ijUj rather than just the linear form as 
done here, study representation functions that are analogs of our symmetric 
Krawtchouk matrices. Their L^^p is our symmetric representation scaled by 
multinomial factors. In addition they suggest further interesting generaliza- 
tions beyond symmetric tensors. 
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